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Abstract 

The WIND Code, a general purpose Navier-Stokes 
solver, has been utilized to obtain supersonic external 
flowfield CFD solutions over an axisymmetric, para- 
bolic forebody with comparisons made to wind tunnel 
experimental results. Various cases have been investi- 
gated at supersonic freestream conditions ranging from 
Mach 2.0 to 3.5, at 0° and 3° angles-of-attack, and with 
either a sharp-nose or blunt-nose forebody configur- 
ation. Both a turbulent (Baldwin-Lomax algebraic tur- 
bulence model) and a laminar model have been 
implemented in the CFD. Obtaining the solutions 
involved utilizing either the parabolized- or full-Navier- 
Stokes analyses supplied in WIND. Comparisons have 
been made with static pressure measurements, with 
boundary-layer rake and flowfield rake pitot pressure 
measurements, and with temperature-sensitive paint 
experimental results. Using WIND'S parabolized Navier- 
Stokes capability, grid sequencing, and the Baldwin- 
Lomax algebraic turbulence model allowed for 
significant reductions in computational time while still 
providing good agreement with experiment. Given that 
CFD and experiment compare well, WIND is found to 
be a good computational platform for solving this type of 
forebody problem, and the grids developed in 
conjunction with it will be used in the future to 
investigate varying freestream conditions not tested 
experimentally. 

Introduction 

A reusable, Rocket Based Combined Cycle 
(RBCC), single-stage-to-orbit launch vehicle concept 
called GTX (formerly known as Trailblazer) 1 is one that 
will use air-breathing propulsion to aid in reducing the 
required propellant mass fraction and, therefore, cost of 
access to space. Figure 1(a) shows a scaled-down 
model of the vehicle concept that was tested in NASA 
Glenn’s 10-by-10 foot Supersonic Wind Tunnel (SWT) 
with Mach numbers ranging from 2.0 to 3.5 and at 
various angles-of-attack to experimentally evaluate the 
external aerodynamics and mass capture character- 
istics. 2 It was desired to model the experimental set-up 
(as far as the forebody is concerned) and compute the 
external flowfield to a reasonably accurate extent with 
CFD. This was done in order to develop a capability 


that could supply any additional external flowfield 
information (which may eventually affect the inlets) not 
supplied by the experiment. As part of the project 
scope, the use of grid sequencing and a parabolized 
Navier-Stokes 3 (PNS) analysis was investigated in 
order to save computational time and yet still obtain 
high fidelity CFD results. A “simultaneous” integration 
between experimental and CFD efforts was made 
possible because of these capabilities. 

Methods And Procedures 

Experiment 

The forebody hardware had the following parabolic 
contour: r=-0.001 16105 x 2 +0.223959x, where r=radius 
and x=axial distance. Either a sharp-nose tip or a 
0.464 in. radius blunt-nose tip was used during each test. 
This parabolic geometry transitioned smoothly into a 
constant radius cylindrical section from x=96.45 in. and 
aft (referenced from the sharp-nose forebody tip). The 
conditions investigated for the sharp- and blunt-nose 
forebodies included Mach 2.0, 2.5, 3.0 and 3.5, with a 
unit Reynolds Number (Re#/ft) of -3.3X10 6 , 2.7X10 6 , 
2.5X10 6 and 2.5X10 6 , respectively. Angles-of-attack (a) 
equal to 0°, ±3°, ±6° and ±9° were investigated at each 
freestream condition, though only 0° and 3° are focused 
on in this paper. All angles are with respect to the 
model’s installation in the 10- by 10-foot SWT. 

During several sharp-nose experimental tests for 
a=0°, a grit was applied at x=3.0 in. in order to trip the 
boundary-layer to turbulent. Sharp-nose forebody 
experimental results with and without application of the 
grit will be presented and compared to the turbulent and 
laminar CFD solutions. 

For the sharp-nose configuration, the CFD results 
can be directly compared with the following 
experimental data: 

1) Static pressure measurements along the forebody, 

2) Boundary-layer rake and flowfield rake meas- 
urements (both measuring pitot pressures) at 
various axial and circumferential forebody 
locations, and 

3) Temperature-sensitive paint (TSP) results. 
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For the blunt-nose forebody, only the boundary-layer 
rake pitot pressure measurements were available and 
could be compared with the CFD solutions. 

Computational Model 

Coordinates for the forebody were supplied and 
incorporated into the gridding process. A rough 
schematic and general description of what was modeled 
in the CFD is shown in Fig. 1(b). The experimental data 
to be used in comparison with the CFD data is relevant 
up to a point just in front of the three propulsion pods’ 
centerbody tips at x=93.33 in. (which are visible in 
Fig. 1(a)). The pods were not modeled in the CFD 
process. However, the grid did model the constant 
cylindrical section up to x=150 in. Therefore, experi- 
mental data obtained circumferentially between the pod 
positions is only valid (in comparisons with the CFD) at 
locations where the first effects from each centerbody’ s 
conical shock would still not be disturbing the forebody 
flowfield. 

In modeling the sharp-nose forebody configuration, 
a 3-D 1 80° single zone, structured grid with two levels 
of sequencing was created (modeling x=-25 to 150 in.) 
using ICEM CFD Mulcad. 4 The grid is shown in 
Fig. 2(a). The grid contains 161 grid points axially in 
the /-direction, 149 radially in the ^-direction, and 
61 circumferentially in the ^-direction. Because of the 
way the sharp-nose grid was assembled, the PNS 
capability (for space marching in the axial direction) 
could be used in the WIND solution process. The wall 
spacing at the body was 0.0002 in. (based on y + 
criteria), with 0.1 in. axial spacing at the nose and 
1.0 in. axial spacing along the forebody surface. The 
sharp-nose grid /= 1 boundary condition (B.C.) was set 
to a frozen B.C., as was the j=j max B.C. An outflow B.C. 
was used for /=/ maY . The j= 1 B.C. consisted of a 
singular axis B.C. ahead of the forebody and a viscous 
wall B.C. along the forebody. For this 180° grid, k= 1 
and A-& max were reflection B.C.’s. 

A separate 3-D 1 80° grid was created for the blunt- 
nose forebody configuration, with a 0.464 in. radius 
nose-tip. The grid created for this blunt-nose forebody 
can be seen in Fig. 2(b). To accurately grid the 
forebody in the nose-tip region, a somewhat “polar 
grid” was created, with the /=1 location corresponding 
to the singular axis boundary condition. The grid wraps 
around the nose while traveling in the /-direction. The 
grid dimensioning corresponds to 141x149x61 (for two 
levels of sequencing). The wall spacing was equal to 
0.0002 in., with 0.06 in. axial spacing at the nose, and 
1.0 in. axial spacing along the body. An FNS analysis 
instead of a PNS analysis was the solution method used 
for the blunt-nose grid because of the singular axis at 
/=1. Depending on processor availability, for certain 
runs the grid could be split into two zones and run using 
WIND’S parallel processor capability. 


The WIND Code’s capabilities and documentation 
are described in detail in Ref. 5. WIND Version 3.0 was 
used in all CFD calculations. The viscous equations 
were solved during each run of WIND using either a 
turbulent (Baldwin-Lomax algebraic model) or a 
laminar model. WIND was run in standard mode, using 
the ideal gas model along with the Sutherland law for 
the transport properties. An adiabatic wall condition 
was utilized. The parameters specified within the 
WIND input file included total freestream conditions, a 
downstream keyword of pressure extrapolate always , 
specification of the turbulence or laminar model, 
sequencing specification, and PNS- or FNS-specific 
commands. WIND was run mainly in single processor 
mode using either R 10000 (225 MHz) or R 12000 
(400MHz) SGI processors, depending on availability. 

The sharp-nose grid solutions involved running the 
PNS capability to convergence. In addition to 
accurately resolving the boundary-layer, the L2 norm of 
the residuals was seen to drop ~3 orders of magnitude 
(except at the nose-tip) to a value usually below 10~\ 
The blunt-nose grid (with FNS capability) was 
investigated for convergence by monitoring the 
integrated pressure and viscous forces, the boundary- 
layer profile, the flowfield, and usually a drop of 
-3 orders of magnitude in the L2 norm of the residuals. 
The primary grid sequencing used corresponded to 
every other grid point from the grid dimensions quoted 
above (both for sharp- and blunt-nose cases). This 
approach was selected based on a grid sequencing/ 
convergence study (per the methods of Ref. 6) 
described in the Appendix. 

Results And Discussion 

Results from the grid convergence study (see 
Tables 1 and 2) showed the following general trends 
with respect to CPU time. The PNS cases (for sharp- 
nose grids) ran 5 to 10 times faster than the FNS cases 
(blunt-nose grid). Using the Medium grids (every other 
grid point) instead of the Fine grids reduced the run 
time by approximately a factor of 25. Use of the PNS 
capability where possible, use of the coarser (Medium 
grid spacing) grids, and use of an algebraic turbulence 
model (instead of a one- or two-equation model) helped 
greatly reduce the computational time, given that most 
of the cases were run on a single processor. This 
allowed for a greater number of comparisons that could 
be made with experimental results. These comparisons 
will now be discussed. 

Sharp-nose Results 

Normalized Surface Static Pressure Distributions. 
Fig. 3(a) displays the normalized static pressure (P stat k/ 
Ptunnei static) plotted versus the axial distance from the 
nose-tip for the Mach 2.5, a=0° test case. Turbulent 
CFD results are plotted with the corresponding experi- 
mental (without grit-applied) static pressure data. The 


N AS A/TM— 200 1 -2 1 0568 


2 


CFD results are seen to compare well with the 
experimental data. At ot=0°, it is expected that all static 
pressures at a given axial station along the circum- 
ferential (0-direction) should be equal (i.e., no pitch or 
yaw). The figure shows that, in general, this is the case 
here. 0 is measured, looking upstream, traveling 
counter-clockwise (or clockwise because of symmetry) 
from the top of the model’s forebody to the bottom (see 
Fig. 1(b)). The static pressure is highest at the nose-tip, 
after which it steadily decreases as the flow expands 
over the parabolic forebody surface. The slight 
waviness in the CFD solution near the nose-tip region 
may be due to the way the PNS solution method 
handles the coarser grid at this “discontinuity.” 

For Mach 2.5 (a=3°), Fig. 3(b) shows CFD and 
experimental results coincide fairly well. When at 
angle-of-attack, higher static pressures are expected on 
the windward side (0=180° line) of the forebody as 
compared to the 1207240° lines, the 907270° lines and 
the leeward side (0° line). The 0° line along the 
forebody surface experiences the lowest pressures. The 
experimental 90° and 270° values, as well as the 120° 
and 240° values again line up fairly well. Note that the 
static pressure distribution for the oc=0° case (Fig. 3(a)) 
would lie approximately centered in between the a=3° 
distributions presented for the 0° to 180° circum- 
ferential range in Fig. 3(b). 

In general, the trends at the other Mach numbers 
investigated are the same as for Mach 2.5, though the 
Mach 3.0 and 3.5 cases both showed several 
discrepancies. As is seen in Fig. 4(a), the Mach 3.0 
solution showed the experimental normalized static 
pressure distribution consistently offset below the CFD 
solution. As is seen in Fig. 4(b), the Mach 3.5 static 
pressure distributions showed good agreement at the 
front of the forebody, but the experimental static 
pressures were found to be lower than the CFD static 
pressure distribution line further downstream. The 
source of the discrepancies could be attributed to 
increases in velocities in the SWT not accounted for in 
the CFD. Axial (in addition to planar) wind tunnel 
gradients on the order of AMach~0.04 spanning over 
the volume occupied by the model were found in the 
Mach 3.0 and 3.5 calibration curves (at high Re#/ft). 
This is important since all CFD cases were run exactly 
at 2.00, 2.50, 3.00 and 3.50. As will be discussed in the 
boundary-layer rake section, even the small Mach 
changes mentioned above could help account for the 
discrepancies found in the results. 

Normalized Flowfield Rake Pitot Pressure Contour 
Plots at x=98.95 in. (0=5° to 35°.175° to 145°) 

A comparison of CFD and experimental flowfield 
rake data obtained at x=98.95 in. is now made. Located 
near a plane slightly aft of the centerbody tips, flowfield 
rakes were used to measure the variation of the pitot 
pressure circumferentially on the leeward and wind- 
ward sides of the forebody. It should be noted that 


experimentally there were only four rakes (0=5°, 15°, 
25° and 35°) used for each run, but the model was run 
at both positive and negative angles-of-attack. This 
allowed the data to be compiled as shown in Fig. 5 
using data from eight rakes (adding data at 0=145°, 
155°, 165° and 175°). Figure 1(b) gives a general idea 
of where these eight rakes (at x=98.95 in.) are located. 
The variable plotted in Fig. 5 is the normalized pitot 
pressure (P pit ot/Ptunnei total)- The CFD solution used in the 
comparison is that obtained assuming a turbulent 
boundary-layer over the entire forebody. The two left- 
hand side contour quadrants (top and bottom) represent 
the experimental data, and they can be directly com- 
pared with the CFD results that have been placed to the 
right-hand side. Refer to the visual aids to the right of 
the quadrant plots to better understand the data being 
presented. The surface of the forebody has been 
depicted in each quadrant figure with a dashed line. 

The particular cases shown in Fig. 5 are the 
Mach 2.5 comparison for a=3° (Fig. 5(a)) and the 
Mach 3.5 comparison for a=3° (Fig. 5(b)). Overall, 
comparing CFD to experiment, the agreement is very 
good for all the conditions tested. The bottom 
(windward) rakes show higher gradients closer to the 
forebody surface as compared to the top (leeward) 
rakes, as may be expected when at an angle-of-attack. 
WIND’S PNS capability is seen to be capturing the 
general features of the flowfield quite well. The trends 
at the other Mach numbers are very similar to those just 
discussed. Generally, the windward side experienced 
slightly higher pitot pressures than the leeward in the 
outer “freestream” region of the external forebody flow. 

Temperature-Sensitive Paint (TSP) Data . 

Figure 6(a) shows experimental TSP data for each 
of the four Mach number cases at a=0° (without grit 
application). The equilibrium surface temperature can 
be inferred by the TSP surface profile distributions. 
This plot can help characterize the transition from a 
laminar to a turbulent boundary-layer occurring along 
the forebody. For Mach 2.0, 2.5 and 3.0, similar 
temperature distributions are found where the 
minimum-maximum-minimum peaks are occurring at 
the same axial locations along the forebody. The 
general location for this “transition” region is thus 
occurring from x~20 to ~35 in. The sudden increase in 
temperature is known to correspond to the beginning of 
transition. 7 Though the Mach 3.5 curve differs from the 
other three, it appears the peak occurs slightly aft of the 
peak for the other Mach numbers. Taken as a whole, 
however, this data shows transition is occurring at 
roughly the same location, which may be expected 
since the Re#/ft is approxi-mately the same for all four 
cases. 

Figure 6(b) shows the same four flow conditions, 
but now the grit has been applied at x=3.0 in. This is 
evident from the “blip” in the temperature distribution 
seen at this location. Comparing with and without grit 
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application, the temperature distributions have changed. 
Both Mach 2.0 and 2.5 (with grit applied) show 
approximately a constant temperature level that is near 
that seen on the aft portion of the forebody in the no grit 
applied case (indicating the grit has tripped the 
boundary-layer to turbulent further upstream). Though 
Mach 3.0 and 3.5 with grit differ, transition to turbulent 
is seen to have occurred further upstream than in the no 
grit cases. This experimental TSP data (along with the 
experimental boundary-layer rake data) will be com- 
pared with the turbulent and laminar CFD results in the 
next section to investigate this transition further. 

Normalized Boundary-Layer Rake Pitot Pressure Profiles 

In the experiment, four boundary-layer rakes were 
placed in various positions axially and circumferentially 
along the forebody surface. For each test condition, the 
model was run at both a positive and negative angle-of- 
attack. Therefore, experimental data can be reduced to 
eight rake locations using symmetry considerations and 
positive/negative angle-of-attack information (see 
Fig. 1(b)). The rake locations through symmetry 
considerations reduce to the following: (x=32.5 in., 
0=58°, 122°), (x=53.94 in., 0=9°, 171°), (x=75.38 in., 
0=87°,93°), and (x=100.33 in., 0=10°, 170°). The 
distance normal to the surface versus the normalized 
pitot pressure (P pil0 t/Ptunnei total) will be plotted. Turbulent 
and laminar sharp-nose grid CFD solutions have been 
obtained for all the cases investigated and are compared 
with experimental results in this section. Plotting both 
CFD solutions together with the experimental data 
allows for a “first order” estimate of the type of 
boundary-layer (laminar, transitional, or turbulent) over 
the entire forebody. Though this is not quite accurate 
(since the downstream boundary-layer is dependent on 
what the boundary-layer characteristic has been 
upstream), the comparison will be shown to be good 
enough for the purposes of this study and provides 
“bounds” on the type of flow. 

Figure 7 shows the four boundary -layer rakes for 
the Mach 2.0 case at oc=0°, with additional experi- 
mental data included for the case where grit was 
applied to the sharp-nose forebody configuration to trip 
the boundary-layer to turbulent. Over the forebody, the 
experimental boundary-layer data without grit applied 
shows good agreement with the turbulent CFD data at 
the three downstream rake locations. However, the 
boundary-layer rake located at x=32.5 in. shows the 
data does not entirely lie on the turbulent solution 
profile, but rather has a slightly thinner profile. This 
agrees with the no grit TSP results presented 
previously, which showed that the turbulent boundary- 
layer was not obtained until somewhere near this rake 
position. 

The assumption made in the flowfield rake section 
discussed previously (in regards to presenting the 
experimental flowfield rake data at x=98.95 in. with 
only the turbulent CFD data) is quite sufficient for 


sharp-nose forebody comparisons. Note that laminar 
and turbulent CFD static pressure distributions were 
very identical. Hence, this was the reason for only 
plotting experimental static pressure data with the 
turbulent CFD data. 

Comparing the Mach 2.0 results with the grit- 
applied experimental data, it can be seen in Figs. 7(a-d) 
that the profiles lie almost on top of the turbulent CFD 
solution profiles. The profiles are seen to be slightly 
thicker than profiles without grit. Transition has 
occurred further upstream than the case without grit. 
Therefore, agreement with the turbulent CFD profile is 
much better for the case with grit. 

It should also be noted, though not presented here, 
that all eight experimental rake profiles for the a=3° 
case (without grit) lie on top of their respective 
turbulent CFD profiles. This is the case for all cases 
tested with the sharp-nose grid. Therefore, there is very 
good experimental and CFD agreement at angle-of- 
attack for the sharp-nose forebody, and WIND (with all 
the assumptions made) has done a good job of captur- 
ing the pertinent flow phenomenon. This suggests that 
transition from laminar to turbulent occurs further 
upstream when at angle-of-attack than is seen for the 
oc=0°. As would be expected at an angle-of-attack, the 
leeward side rakes, especially those closest to the top of 
the forebody (x=53.94 and 100.33 in.), display a thicker 
boundary-layer profile than the windward side rakes at 
the same axial locations. In addition, as the flow 
proceeds downstream over the forebody, the effects of 
expansion can be seen to be taking place, with 
“freestream” pitot pressures continually decreasing. 
Also note that the pitot pressure does not reach a 
constant value outside the “boundary-layer” in these 
plots because of the “conical” nature of the flow. 

Mach 2.5 and 3.0 sharp-nose forebody experi- 
mental data showed very similar trends to the Mach 2.0, 
a=0° and 3° results just discussed (i.e., experiment and 
turbulent CFD pitot pressure profiles don’t line up 
exactly until after x=32.5 in. for a=0°, and they line up 
at all locations measured for a=3°). However, Mach 3.5 
actually showed that the experimental boundary-layer 
data without grit does not line up with the turbulent 
CFD profile until somewhere after x=53.94 in. Transi- 
tion appears to have occurred further back as compared 
to the other Mach number cases (as might be 
ascertained from the TSP data in Fig. 6(a)). 

Mach 3.5 pitot pressures in the “freestream” (for 
the experiment and CFD) were offset from one another 
at various axial stations. Higher “freestream” pitot 
pressures were obtained for the experimental data as 
compared to the CFD data at upstream locations, while 
lower “freestream” experi-mental pitot pressures were 
obtained at the downstream locations. The trend agreed 
with 10- by 10-foot SWT data showing an increasing 
freestream Mach number with axial travel downstream 
in the tunnel for certain conditions (and therefore a 
decreasing “freestream” pitot pressure). Several Mach 
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number perturbation CFD runs (on the order AMach 
-0.04) showed the discrep-ancies between experimental 
and CFD boundary-layer rake “freestream” values 
could be explained by these slight Mach number 
gradients. 

Blunt-nose Results 

Numerous FNS runs of the WIND code have been 
completed for the blunt-nose Medium grid (every other 
grid point) at a=0° and at a=3° (M = 2.0, 2.5, 3.0, and 
3.5). A summary of the results for these cases follows. 

Comparing boundary-layer profiles, trends at Mach 
2.0 and Mach 2.5 (for a=0°) are somewhat similar to 
the sharp-nose Mach 2.0, 2.5 and 3.0 (a=0°) forebody 
trends. Figures 8(a-d) shows the Mach 2.5, a=0°, blunt- 
nose boundary-layer rake profile comparisons (all 
without grit) over the forebody length. At x=32.5 in., 
the experimental data actually is seen to lie closer to a 
laminar boundary-layer profile than a turbulent one. 
However, further downstream at x=53.94 in., the 
experimental profiles are lying on top of the turbulent 
CFD solution profiles. Generally speaking for these 
cases just mentioned, the thickness for both the sharp- 
and blunt-nose profiles, as well as the normalized 
“freestream” pitot pressure values, agree quite well with 
each other (station-to-station) over the forebody 
surface. 

Mach 3.0 and 3.5 boundary-layer rake results, 
however, show that the experimental data (no grit was 
applied in any of the blunt-nose cases) do not lie 
directly on the turbulent CFD solution until further 
downstream than had been seen in the sharp-nose 
forebody solutions. For the Mach 3.5, a=0° solution, 
the experimental profile lies atop the laminar CFD 
profile even up to x=53.94 in. Therefore, transition is 
occurring downstream of this station for the blunt-nose 
condition, whereas this type of definitive statement 
could not be made for the corresponding sharp-nose 
cases (where the upstream experiment profiles were 
found only slightly offset from the turbulent CFD 
solutions). 

For a=3°, the trends are actually slightly different 
for several of the blunt-nose cases as compared to the 
sharp-nose cases. As before for the sharp-nose cases at 
Mach 2.0 and 2.5, the experimental profiles lay nearly 
atop the turbulent CFD profiles at all the rake locations. 
However, the Mach 3.0 and 3.5 experimental profiles 
(up to and including x=53.94 in. on the windward side) 
lie either somewhere in between turbulent and laminar 
CFD profiles or directly on top of the laminar profile 
(see Fig. 9 for Mach 3.0). This was not seen in the 
sharp-nose cases at any Mach number when the 
forebody was at any angle-of-attack. 

In summary, it may be deduced that, in general, 
transition has occurred at axial locations ahead of where 
the experimental profile finally matches the turbulent 
solution. When taking all boundary-layer rakes as a 
whole, this overall view of the data can be compared to 


the experimental TSP data (as was shown in Fig. 6). 
The conclusions drawn from the CFD-combined experi- 
mental normalized pitot pressure plots show good 
general agreement with the TSP data conclusions, and 
provide some general idea of boundary-layer transition 
location and characteristics. Overall, the PNS capability 
is seen to have done a good job in returning worthwhile 
results that agree well with experiment. 

Summary And Conclusion 

In comparisons involving NASA Glenn’s 10- by 
10-foot Supersonic Wind Tunnel experimental results, 
very good agreement between WIND Code CFD and 
experiment has been found. Using WIND’S PNS cap- 
ability where possible, the coarser grids, and the 
Baldwin-Lomax algebraic turbulence model allowed 
for great reductions in the computational time while 
still maintaining good agreement with the experimental 
results. Running CFD concurrently with experimental 
testing allowed for one’s information to feed directly 
into the planning for the other. Future use of the CFD 
models to obtain reasonably accurate information that 
was not obtained with the experiment is now an option. 
The WIND code has been found to be a good solver for 
the applications set forth in this study. 

Appendix 

Grid Convergence Study 

A grid convergence study for the sharp-nose and 
blunt-nose grids has been performed based on the 
procedures outlined in Ref. 6. The total drag coefficient 
(pressure plus viscous drag coefficient) was the 
parameter used for comparison in these studies. 

First, for the sharp-nose grid, the Mach 2.5 (a=0°) 
flowfield was computed on three grids, each with twice 
the number of grid points in the /-, y-, and k- directions 
as the previous grid (i.e., Coarse, Medium and Fine grid 
spacing). The Fine grid was chosen as the reference 
grid and is referred to as grid 1, the Medium grid is 
referred to as grid 2, and the Coarse as grid 3. 

The Grid Convergence Index (GCI), similar to the 
Richardson error estimator and based on the 
generalized theory of Richardson Extrapolation, 
provides a calculated error band on the solution 
functionals (total drag coefficient in this case). It also 
provides a measure of uncertainty for grid convergence. 
Though two levels of sequencing (three grids) were 
utilized with a constant grid refinement ratio of 2, the 
GCI was calculated using a Factor of safety (F s ) equal 
to 3 in order to be conservative. Table 1 presents the 
sharp-nose grid’s relevant information for the grid 
convergence study. Defining “/?” (the order of 
convergence) per Ref. 6, it is calculated by using the 
solution functionals obtained from the three grids and 
the grid refinement ratio. For the sharp-nose grid, using 
Table 1 values, p was calculated as -1.74. This shows 
the solution is not quite at a second-order level of 
convergence. This may be expected since, with the PNS 
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capability, WIND automatically changes to a first-order 
explicit operator in the /-direction (marching direction) 
only, while keeping second-order operators in the other 
directions. The Fine grid GCI (GCI 12 ) was calculated as 
0.05 percent, while the Medium grid GCI (GCI 23 ) was 
calculated as 0.17 percent. The solutions were found to 
be in the asymptotic range [GCI 23 /(GCI 12 *2 P ) » 1.00]. 
These values were deemed good enough for the 
purposes of this study. The Medium grid has been used 
throughout this paper for the sharp-nose forebody cases. 

For the blunt-nose grid, the Mach 2.5 (a=0°) 
flowfield was also computed on three grids, each with 
twice the number of grid points in the /-, 7 -, and 
A-directions as the previous grid (i.e., Coarse, Medium 
and Fine grid spacing). Using the same procedure as 
described above, Table 2 presents the relevant blunt- 
nose grid information for the grid convergence study. 
The order of convergence p was calculated as -1.52. 
This shows that the solution is not at a second-order 
level of convergence, where FNS calculations for the 
blunt-nose case use second-order explicit operators in 
all directions. The difference between theoretical and 
calculated order of convergence may be due to grid 
stretching, grid quality, nonlinearities in the solutions 
(such as the presence of shocks) and possibly other 
factors . 6 The Fine grid GCI (GCI 12 ) was calculated as 
0.38 percent, while the Medium grid GCI (GCI 23 ) was 
calculated as 1 .09 percent. The solutions were found to 
be in the asymptotic range [(GCI 23 /(GCI, 2 *2 P )) ~ 1.00]. 
These values were deemed good enough for the 
purposes of this study. The Medium grid has been used 
throughout this paper for the blunt-nose forebody cases. 
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Table 1 . — Sharp-nose Forebody Grid Convergence Study Results - 

Mach 2.5, <x=0° 

Grid Number 

Grid Spacing 
(Normalized by 
finest grid spacing) 

Total CPU time 
(Corrected for 
machine use) 

Total Drag 
Coefficient 

Calculated order 
of convergence 
(P) 

Fine ( 1 ) - 000 sequencing in WIND 

1 

-75,500 sec 

0.108639 


Medium (2) - 1 1 1 sequencing in WIND 

2 

-3,000 sec 

0.108596 

1.74 

Coarse (3) - 222 sequencing in WIND 

4 

-200 sec 

0.108452 


Table 2. — Blunt-nose Forebody Grid Convergence Study Results - 

Mach 2.5, a=0°. 

Grid Number 

Grid Spacing 
(Normalized by 
finest grid spacing) 

Total CPU time 
(Corrected for 
machine use) 

Total Drag 
Coefficient 

Calculated order 
of convergence 
ip) 

Fine ( 1) - 000 sequencing in WIND 

1 

-446,500 sec 

0.107664 


Medium (2) - 1 1 1 sequencing in WIND 

2 

-17,000 sec 

0.107411 

1.52 

Coarse (3) - 222 sequencing in WIND 

4 

- 2,000 sec 

0.106686 
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x=98.95 in 



(a) (b) 

Figure 1 - (a) Scaled-down GTX Model, (b) Schematic and general description of forebody used for CFD. 




(a) (b) 

Figure 2 - Isometric views, (a) Sharp-nose forebody grid, (b) Blunt-nose forebody grid. 
(Every grid point for the finest grid is shown.) 
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P(statlc)/P(tunnel static) P(statlc)/P(tunnel static) 




X - (distance from nose tip) x - (distance from nose tip) 


(a) (b) 

Figure 3 - Sharp-nose Forebody Normalized Static Pressure Distribution, (a) Mach 2.5, oc=0°, 

(b) Mach 2.5, a=3°. 



10 20 30 40 50 60 70 80 

X • (distance from nose tip) 



X - (distance from nose tip) 


(a) (b) 

Figure 4 - Sharp-nose Forebody Normalized Static Pressure Distribution, (a) Mach 3.0, a=0°, 

(b) Mach 3.5, a=0°. 
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Distance from vehicle centerline (inch) 

(a) 


Experiment 


CFD 





Leeward 


Windward 


Symmetry 

Plane 



Distance from vehicle centerline (inch) 

(b) 


Figure 5 - Sharp-nose Flowfield Rake Normalized Pitot Pressure Contour Plots for a=3° at x=98.95 inches 

(experiment vs. turbulent CFD), 

(a) Mach 2.5, a=3°, 

(b) Mach 3.5, a=3°. 


NASA/TM — 200 1 -2 1 0568 


9 


Normalized Temperature 




(a) (b) 

Figure 6 - Sharp-nose normalized temperature distribution over part of forebody for ot=0° cases, 
(a) No grit applied, (b) With grit applied at x=3.0 inches. 



Figure 7 (a-d) - Sharp-nose Boundary-Layer Rake Normalized Pitot Pressure Profiles along the forebody at various 
axial stations (turbulent and laminar CFD vs. experimental (with and without grit)), 

Mach 2.0, a=0°, 

(a) x=32.5 inches, (b) x=53.94 inches, (c) x=75.38 inches, (d) x=100.33 inches. 
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Figure 8 (a-d) - Blunt-nose Boundary-Layer Rake Normalized Pitot Pressure Profiles along the forebody at various 
axial stations (turbulent and laminar CFD vs. experimental (all without grit)), Mach 2.5, a=0°, 

(a) x=32.5 inches, (b) x=53.94 inches, (c) x=75.38 inches, (d) x= 100.33 inches. 



^pito/^tunnel total 


Figure 9 - Blunt-nose Boundary-Layer Rake Normalized Pitot Pressure Profiles along the forebody at various axial 

stations, Mach 3.0, a=3°. 
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